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A transport approach including a leakage effect for J ftp's 
in the transverse phase space is used to calculate the ratios 
between the J/ip transverse momentum distributions in sev- 
eral centrality bins for Pb-Pb collisions at E/A = 160 GeV. 
From the comparison with the CERN-SPS data, where the 
centrality is characterized by the transverse energy Et, the 
leakage effect is extremely important in the region of high 
transverse momentum and high transverse energy, and both 
the threshold and the comover models can describe the ra- 
tio well for all centrality bins except the most central one 
(E t < 100 GeV), for which the comover model calculation is 
considerably better than the threshold one. 
PACS : 25.75. - q, 12.38. Mh 

Relativistic heavy ion collisions may be the only way to 
create the extreme condition for producing a new state 
of matter — Quark-Gluon Plasma (QGP). Among the 
possible signatures suggested to identify the existence of 
QGP, the J/ip suppression proposed by Matsui and Satz 
[1] is considered as a direct one of deconfinement phase 
transition. The discovery of anomalous J/ip suppression 
[2] in Pb - Pb collisions in 1996 has been one of the high- 
lights of the research with relativistic nuclear collisions 
at SPS. Since different models [3,4], with and without 
the assumption of QGP, can describe the observed sup- 
pression, it is still not yet clear if the suppression means 
the discovery of QGP. Recently, the NA50 collaboration 
has presented the final analysis of the J/ip transverse 
momentum distribution [5] in Pb - Pb collisions at E/A 
= 160 GeV. The mean squared transverse momentum 
(pf) as a function of the centrality of the collisions first 
increases and then becomes flat when the centrality in- 
creases. The understanding of this saturation of (pf) in 
central collisions needs the consideration of time struc- 
ture [6,7] of the anomalous suppression. With the idea 
of leakage, which has already been considered more than 
10 years ago [1,8,9] (c.f. also more recent works [10,11]), 
the anomalous suppression is not an instantaneous pro- 
cess, but takes a certain time. During this time the J/ip's 
with high transverse momenta may leak out of the source 
of the anomalous suppression. As a consequence, low p t 
J/ip's are absorbed preferentially, and the mean pt of the 
survived J/ip's in central collisions will not drop down 
but become flat. From the comparison of the theoretical 
calculations [7] using two different models, the thresh- 
old [12,13] and comover [14,15] models, with the data, 
the average anomalous suppression time tA is extracted 
to be 3 — 4 fm/c for central Pb - Pb collisions. Obvi- 
ously, the leakage mechanism affects mainly the high p t 



J/ip's, the mean values (p t ) < 1.2 GeV/c and (pf) < 2 
(GeV/c) 2 for Pb - Pb collisions can not reflect the leak- 
age effect sufficiently. In order to see the leakage effect 
more clearly, we calculate in this paper the J ftp pt dis- 
tributions in different centrality bins, and compare them 
with the data of NA50. We will focus our attention on 
the region 1 < p t < 3 GeV/c for central collisions where 
the leakage effect and anomalous suppression are both 
important, and the error bars of the data are small. 

J/ip's produced in a nucleus-nucleus collision experi- 
ence first normal suppression via inelastic J ftp — N colli- 
sions and show normal broadening of (pf) above (pf)jvjv 
via gluon rescattering [4]. A gluon from a projectile nu- 
cleon and a gluon from a target nucleon fuse to form a 
J/ip at a point with space coordinates (s, za) in the rest 
frame of the target A and (6 — s,zb) in the rest frame 
of the projectile B, where b is the impact parameter for 
the BA collision. On its way out, the J/tp experiences 
the thickness Ta(s, za, oo) and Tg(b — s, — oo, zb) in the 
nuclei A and B, and is suppressed with an effective ab- 
sorption cross section a a b s , where the thickness function 
T is defined as T(s, Z\, 22) = J z ^ dzp(s, z) with pa,b be- 
ing the density of nuclei A and B. Before the fusion, 
the two gluons traverse thicknesses Ta{s, —00, z a) and 
Tb(6 — s, zb, 00) of nuclear matter in A and B, and ob- 
tain additional transverse momentum via g N collisions. 
Neglecting effects of formation time [16], one has after 
the normal suppression the J / tp distribution function /jv 
in the transverse phase space at given 6, 



(1) 



fN{pt\b, s) = / dz A dz B p A (s, z A )pB{b - s, z B ) x 



,—Oabs (T a(s, z A ,oo)+T B {b-s,- 00, Z B )) _ 



(P 2 t) 



-Pt/(PI)N 



N 



with 

(pt) N (b,s,z A ,z B ) 



{Pt )nn 



dgNPo 1 [Ta(s, -00, z a ) + T B (b~ s, z B ,oo) 



(2) 



where we have assumed a Gaussian p t dependence. The 
constants <j a bs and a g pf are usually adjusted to the data 
from pA collisions where J/ip's experience only normal 
suppression. 

Anomalous suppression is attributed to the action on 
the J/ip by the mostly baryon free phase of partons 
and/or hadrons which is formed after the nuclear overlap. 
We define t = as the time when the normal suppression 
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and normal generation of (p 2 ) have ceased, and describe 
the time evolution of the J/ip density f(p t ,t\b 7 s) in the 
transverse phase space by a transport equation [7] 



dt 



+ vt-V 8 f = -af 



with the initial condition 

f(p u 0\b,s) = f N (p t ,b,s) . 



(3) 



(4) 



The second term on the left handside of the transport 
equation arises from the free streaming of the J/ip with 

transverse velocity v t — p t / ^jp\ + m 2 ^, and the loss 

term is due to anomalous suppression. The function 
a{pt,t\b, s) contains all details about the mechanism of 
anomalous suppression. 

The transport equation (3) can be solved analytically 
with the result [7] 

f{p u t\b,s) =e-f^ a &Sfa-W-* )) f N &\$,?-€ t t) . 

(5) 

If the anomalous suppression ceases at time tf, a(pt,t > 
tf\b,s) = 0, the observed yield Y as a function of trans- 
verse momentum p t and transverse energy E t can be re- 
lated to the phase space density / via 

Y(pt,E t ) = 

J d 2 bd 2 sd 2 p t P(E t \b)f(p t ,t f \b, s)5(Pt - ^Pl+P 2 y) ■ (6) 

Here, P(E t \b) describes the distribution of transverse en- 
ergy E t in events with a given impact parameter b. It is 
chosen as a Gaussian distribution [12] 



P(Et\b) 



g 2q^aN p (b) (7) 



yj2irq 2 aN p (b) 
with the number of participant nucleons defined as 

N p (b) = J d 2 sn p {b,s) , 

n p (b, s) = T A (s, -oo, oo) (l - e -°™T B (i-3,- 00 , 00 )^ 

+ T B (b- s, -oo.oo) (l - e-"™ 7 "' 8 '- 00 ' 00 ' 
(E t )(b) = qN p (b) , 



(8) 



and a = 1.27, q = 0.274 GeV and a NN = 32 mb [12] for 
Pb - Pb collisions at SPS energy. 

In order to see clearly the leakage effect which is im- 
portant for the J ftp's with high transverse momenta pro- 
duced in the events with high transverse energy, we con- 
sider the ratio between the J/ip p t distributions in the 
E t bin i and in the E t bin j, 



Ri/j{Pt) 



f AEl dE t Y( Pu E t )/N^ 
J AEi dEtY(p u Et)/N^ 



(9) 



To compare the ratio with the experimental data [5] , we 
have normalized the distributions to the number of Drell- 
Yan pairs in the same E t bins, 

N AE' = I dE t I d 2 bd 2 sdz A dz B P(E t \b) x 

* MEj J 

PA(s,z A )p B (b-s,z B ). (10) 

Since the physical origin of anomalous suppression is 
not yet clear, we calculate the ratio R for two rather 
different scenarios: Instantaneous absorption involving a 
threshold in the energy density and continuous absorp- 
tion by comovers. 

We begin with the threshold model and use the ap- 
proach by Blaizot et. al. [12] for simplicity. The suppres- 
sion function in this approach can be recovered within 
the transport formalism (5) by setting 



a (pi, t\b, s) = a 9(n p (b, s) - n c )5{t) 



(11) 



and taking the limit a — > oo. J ftp's are totally and in- 
stantaneously destroyed if the energy density which is di- 
rectly proportional to the participant density n p is larger 
than a critical density, and nothing happens elsewhere. 
While the threshold approach successfully describes the 
data in the full E t range of the J/ ip suppression after the 
only one free parameter, n c , is adjusted, the predictions 
for (p 2 )(E t ) are significantly below the data, especially 
in the high E t region. 

It is difficult to understand that the anomalous sup- 
pression happens instantaneously at time t — 0. The 
sudden suppression (11) was changed to a continuous 
suppression between times to and t\ [7], 



a(p t ,t\b,s) 



a 



t\ — t a 



6(n p {b,s) -n c )6>(ti -t)6(t-t ) 



(12) 



From the comparison [7] of the yield Y(E t ) and mean 
transverse momentum (p 2 )(E t ) with the data, the aver- 
age anomalous suppression time t A was found to be a 
function of E t . It increases from to about 4 fm/c when 
E t increases from to the maximum (~ 140 GeV) for Pb 
- Pb collisions. 

With the parameters a a b s = 6.4 mb and n c = 
3.75 fm -2 [12] for the suppression and (p1)nn = 1-11 
(GeV/c) 2 and a gN = 0.081 (GeV/c) 2 fm" 1 [5] for the 
mean transverse momentum, we calculated the ratios 
R-i/jiPt)- We also accounted for the transverse energy 
fluctuations [12] which have been shown to be signifi- 
cant for the explanation of the sharp drop of J/ip sup- 
pression in the domain of very large E t values, by re- 
placing n p by jj^yfip where (E t ) is the mean transverse 
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energy at given b. For the comparison with the data 
[5] of Pb - Pb collisions, we consider five E t bins with 
AEl = [5, 34], [34,60], [60,80], [80,100], [100,140] GeV 
corresponding to i = 1, 2, 3, 4, 5, respectively. The 
theoretical results Ri/i(pt) with i = 2, 3, 4, 5 together 
with the data are shown in Fig. 1. The dashed lines in- 
dicate the results with normal suppression only. Since 
the nuclear absorption deviates from the data and be- 
comes saturated above E t ~ 40 GeV [3,4], the calculated 
ratios Ri/i are always above the data for any values of 
p t and in any E t bin. The dotted lines are calculated 
with the original Blaizot approach (11) without consid- 
ering leakage effect (v t = in (5)), and the solid lines are 
obtained by taking the average anomalous suppression 
time tA = 0, 1, 2, 3, 4 fm corresponding to the E t bins 
i = 1, 2, 3, 4, 5, respectively. It is clear that the leak- 
age effect becomes more and more important when the 
anomalous suppression time increases and/or the trans- 
verse momentum p t increases. While the instantaneous 
threshold approach, namely the dotted lines, can describe 
the data in the region p t < 3 GeV/c when E t is not high 
enough, the calculated ratio R5/1 deviates from the data 
clearly. Due to the leakage effect reflected in the trans- 
port approach (5), high transverse momentum J/ip's may 
escape the anomalous suppression. This is the reason 
why the solid line is always above the dotted line in any 
E t bin. When the degree of centrality is low, the system 
is difficult to reach the critical density n c , the anoma- 
lous suppression is not important, and then the leakage 
effect is weak. With increasing centrality, more anoma- 
lous suppression happens and more J/^'s leak out of the 
suppression region. For the E t bins 1 and 2, almost no 
leakage happens, the solid and dotted lines coincide for 
i? 2 /i- For -R3/1 the leakage is still not important. How- 
ever, when the centrality increases further, the leakage 
effect leads to a remarkable difference between the solid 
and dotted lines. For R4/1 and R5/1 the leakage im- 
proves considerably the threshold approach. Due to the 
big error bars of the data above p t around 3 GeV/c, it 
is difficult to judge the leakage effect in the very high p t 
region from the comparison with the data. 

We now turn to the comover model. The comoving 
partons and/or hadrons lead to a continuous J ftp absorp- 
tion of long duration due to inelastic collisions with the 
comoving particles. As a representative model, we use 
the approaches by Capella et al. [14,15] and Kharzeev et. 
al. [17]. The absorption term a in the transport equation 
(3) takes the form [7] 

r- +iiT -> n c (b,s) fn c (b,s) \ 
a(p t ,t\b, s) = cr co 6 1 -t]6(t- t ) , 

t \ n f J 

(13) 

where a co is the cross section of J/i/>-comover interaction, 
n c the comover density at initial time to. The absorption 
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FIG. 1. The ratio between the J/tp pt distributions in the 
E t bin i (i — 2, 3, 4, 5) and in the first E t bin (i = 1). 
The distributions are normalized to the Drell-Yan number in 
the same Et bins. The five Et bins correspond to the Et 
range AE t = [5,34], [34,60], [60,80], [80,100], [100,140]. 
The data are from NA50. Dashed, dotted and solid lines 
are calculated within the threshold model. The dashed lines 
show the result of normal suppression alone, the dotted lines 
include the anomalous suppression but without leakage effect, 
and the solid lines include also the leakage effect with the 
average anomalous suppression time tA = 0, 1, 2, 3, 4 fm/c 
corresponding to the E t bin i = 1, 2, 3, 4, 5, respectively. 

by comovers starts at to and ends at t\ when the co- 
mover density n c (b, s)to/t\ has reached a value n/ inde- 
pendent of b and s. Unlike the threshold scenario for 
which we have to introduce a time structure (12) to dis- 
cuss the leakage effect, the comover scenario contains a 
definite time structure for anomalous suppression and it 
is not necessary to introduce any new parameter. Tak- 
ing into account the transverse energy fluctuations [14] 
by replacing n c by j^n c , and the transverse energy 
loss [15] induced by the J/tjj trigger by rescaling (Et) 
by a factor ^ — -, which have been shown to be signif- 
icant for the explanation of the sharp decrease of J/tp 
suppression at E t > 100 GeV, and using the parame- 
ters n co (b,s) = 1.5 n p (b,s),t = 1 fm, n/ = 1 fm~ 2 , 
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c a fcs = 4.5 mb and a co = 1 mb [14,17], the calculated 
(Pt)(E t ) [7] in the transport approach including leakage 
effect agrees well with the data, and the average time 
for comover action is tA — 3.5 fm/c. With the same 
parameters, the calculated ratios Ri/i for i = 2,3,4,5 
are shown in Fig. 2. Again the result with only normal 
suppression deviates from the data. Compared with the 
threshold scenario, the prediction with comover absorp- 
tion but without considering leakage effect has a better 
agreement with the data. Since the comover absorption 
is a continuous process, and it can happen in the events 
with low centrality, the leakage effect starts to change the 
suppression ratio at low E t and is almost E t independent. 
It is clear that the leakage effect leads to a better agree- 
ment between the comover model and the data. From the 
comparison with the threshold scenario, the result of the 
comover scenario looks more close to the experimental 
data. 
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FIG. 2. The ratio between the J ftp pt distributions in the 
E t bin i (i = 2, 3, 4, 5) and in the first E t bin (i = 1). 
The distributions are normalized to the Drell-Yan number in 
the same Et bin. The five Et bins correspond to the Et range 
AE t = [5,34], [34,60], [60,80], [80,100], [100,140]. The data 
are from NA50. Dashed lines (with only normal suppression), 
dotted lines (with anomalous suppression but without leak- 
age effect) and solid lines (with anomalous suppression and 
leakage effect) are calculated within the comover model. 



With the general transport equation proposed in ref. 
[7] which describes the time evolution of the J/^'s during 
anomalous suppression including leakage effect, we calcu- 
lated within the threshold and comover models the ratio 
Ri/i{Pt) between the J/ip p t distribution in the higher 
E t bin i and the same distribution in the lowest E t bin. 
Our motivation is to study the leakage effect in high p t 
and high Et region where the anomalous suppression is 
most important. We found that in this region the leak- 
age effect leads to a better agreement with the data for 
the two models. The results support the idea [7] that the 
anomalous suppression happens on average at 3 — 4 fm/c 
after the normal suppression ends. From the comparison 
of the two very different models, our calculation with the 
mechanism of continuous comover action over a compa- 
rably long time can better describe the ratios for all E t 
bins. 
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